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Abstract 

A high density diquark phase seems to be a generic feature of QCD. If so it 
should also be reproduced by random matrix models. We discuss a specific one 
in which the random matrix elements of the Dirac operator are supplemented 
by a finite chemical potential and by non-random elements which model the 
formation of instanton-anti-instanton molecules. Comparing our results to 
those found in a previous investigation by Vanderheyden and Jackson we find 
additional support for our starting assumption, namely that the existence of 
a high density diquark phase is common to all QCD-like models. 
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I. INTRODUCTION 

The understanding of the QCD phase diagram at finite density and temperature is of 
fundamental interest to investigate phenomena related to heavy-ions collisions and to the 
physics of neutron stars. While first principles lattice calculations are able to take into 
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account the effects of finite temperature, the addition of a baryonic chemical potential to 
the lattice action makes the fermionic determinant complex, which prevents reliable Monte- 
Carlo calculations. Resorting to specific QCD-models or to an asymptotic QCD expansion 
for very large temperature or chemical potential fi is thus necessary to investigate the finite 
density domain. These models allow to study the symmetries of the different regions of 
the phase diagram; of particular interest is the study of the chiral restoration at finite tem- 
perature and/or density. Recently a lot of attention has been focused on so-called "colour 
superconductivity": at large density and low temperature, an arbitrarily weak attractive 
force makes the Fermi sea of quarks unstable with respect to diquark formation and in- 
duces Cooper pairing (diquark condensation) At still higher /x this turns into still 
another, so-called colour- flavour- locked, phase 0. As already stated there exists no really 
reliable QCD technique to analyse these phases, except for completely unrealistically large 
fi. However, all the different approaches used so far yield the same general picture which 
gives credence to the assumption that it is generic. Some of the different approaches fol- 
lowed are N JL-like models ||^ , which can be closely related to models based on the instanton 
phenomenology the Dyson-Schwinger approach 0, and, at very large /i, weak coupling 
expansions These calculations agree in their main conclusions and suggest a fairly large 
gap (on the order of 100 MeV), which could have some observational consequences in neu- 
tron stars and possibly even heavy-ion collisions (for recent reviews see One should 
note however that colour-superconductivity in QCD differs from usual superconductivity 
in the following sense: for large chemical potential, the Renormalisation Group Approach 
leads to an asymptotic dependence of the gap as a function of the QCD coupling constant 
of the kind A ~ exp {—c/ g) [1^. The naive expectation from BCS theory would be in 
contrast A ~ exp {—c/g^). This difference in behaviour is due to the long-range nature of 
the magnetic interaction in QCD: while non-static magnetic modes of the gluon propagator 
are dynamically screened due to Landau damping, static magnetic modes are not screened. 
It leads to a considerable enhancement of the QCD superconducting gap. 

To test the hypotheses that these new phases are indeed a generic feature of QCD and 
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QCD-like models, such that the insufficiencies of the individual approaches are unimportant 
for the existence of these phases, Random Matrix Theory (RMT) is the ideal tool. The 
main idea of RMT is in fact to isolate generic features. In the past RMT was succesfuUy 
used to analyse the eigenvalues of the QCD Dirac operator as obtained in lattice gauge 
calculations. For a recent review see . The perfect agreement between RMT predictions 



and lattice-QCD results led to the general conviction that RMT does apply to QCD, though 
a rigorous formal proof is still missing. 

RMT is especially suited to investigate the chiral phase transition. The Banks-Casher 
relation: 

<qq>= -^p{>^)\x=o (1) 

relates the chiral condensate < gg > to the density p of zero-eigenvalues of the Dirac equa- 
tion and the latter was shown to display universal properties described by RMT. Universal 
spectral correlations of the QCD Dirac operator can also be computed at finite temperature. 
In addition to giving exact results for correlations of eigenvalues, RMT can also be used as 
a schematic model to investigate non-universal quantities. 

Recently Vanderheyden and Jackson have investigated the phase diagram of a QCD-like 
theory with generic 4-quark couplings in the framework of a RMT model [|l2[|l3l with two 
quark flavours. Using a saddle point approximation they derived the thermodynamic po- 
tential for the quark and diquark condensates and studied the competition between chiral 
restoration and diquark condensation as a function of temperature and density. They anal- 
ysed the phase diagram for different values of the coupling constants in the diquark (< qq >) 
and chiral (< qq >) channels and showed that the phase diagram can realize a total of six 
general structures. 

To model the effects of temperature, Vanderheyden and Jackson include the lowest Mat- 



subara frequency only. In a more recent work [Q, they have also investigated the phase 
diagram of A^c = 2 QCD with all Matsubara frequencies included. While the inclusion of all 
Matsubara frequencies doesn't affect the topology of the phase diagram, it eliminates some 



unphysical properties like negative baryonic densities at small /i and the variation of the chi- 
ral condensate with fi. The combination of RMT and the Matsubara formalism is, however, 
only well justified if there are no relevant effects beyond the boundary condition in euclidian 
time. QCD just above the phase transition shows, however, strong correlations between the 
quarks and gluons of a type which can hardly be described by such a simple model. In 
this contribution we therefore try to improve the pioneering paper by Vanderheyden and 
Jackson in this respect by allowing for a more general type of non-random matrix elements. 
We assume that the properties of the lowest Dirac eigenstates are primarily determined by 
instanton and anti-instanton field configurations. In the instanton picture [|15|, spontaneous 



chiral symmetry breaking is assumed to be generated by randomly distributed uncorrelated 
instantons/anti-instantons, which allows for a delocalization of the associated quark quasi- 
zero-modes. The restoration of chiral symmetry at high temperature (or high density) can 
then be realized if the instanton liquid changes from a random ensemble of instantons and 
anti-instantons to a correlated system where instanton and anti-instanton pair to form so- 
called "molecules" (ITH], the precise definition of which is ambiguous but also irrelevant for 



our purposes. The formation of such clusters has been observed on the lattice [|T7| and in 
numerical simulations of the instanton liquid (for more references on this subject, see 
also [§]). Of course, other scenarii for chiral symmetry restoration have been proposed; this 
phase transition can for example also be seen as a Mott- Anderson like transition to an "insu- 



lator state" 191. In this work however, we assume that the dominant mechanism for chiral 



symmetry restoration involves instanton molecular correlations. In ||2^ such a molecular 
model was used to study the chiral phase transition in the framework of RMT at zero den- 
sity. The instanton configuration (and therefore the temperature) was characterized by two 
parameters. It was found that to reduce the value of the chiral condensate by more than a 
factor two, about 95 percent of the instantons and anti-instantons have to pair in molecules, 
which confirmed the decisive role played by molecules formation in the restoration of chiral 
symmetry at finite temperature. In this paper we want to generalize the results of to 



finite chemical potential. That will allow us to investigate the stability of the results found 



in 



1^,0 when the temperature effects are not modelled in the most elementary manner by 
the first Matsubara frequency. Within RMT the properties of the chiral phase transition 
are determined by the interplay between fluctuations (described by the random matrix ele- 
ments) and constant terms (for fixed T and fi) in the Dirac operator. In [0,0 the latter 
were assumed to be the same for all states, while we allow for the possibility that they only 
contribute for a certain fraction of them. We find that for small fi the existence of a phase 
transition is extremely sensitive to this fraction (and thus to the detailed instanton-anti- 
instanton dynamics) while it becomes nearly independent of it for large /i, suggesting that 
the transition to a diquark condensate is indeed a model-independent feature. 



II. FORMULATION OF THE MODEL 



Apart for the parametrization of temperature we follow closely the approach adopted 
by Vanderheyden and Jackson ||12|jl3[| and refer to their original papers for details of the 
derivation of the thermodynamic potential. Our starting partition function (for two quark 
flavours tpi and ^/'2) is the following: 
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i.e., it is the same as in Eq. (1) of |jT3[ but with the substitution ttTItv — ^ D with D being a 
N X N diagonal matrix, for which a fraction a of its diagonal elements are equal to a fixed 
value d and the rest is zero. This form is motivated by the observation that the correlation 
of instantons and anti-instantons (i.e. the formation of molecules) generates such diagonal 
terms. The size of these terms depends on the anti-instanton-instanton separation and d has 
to be interpreted as an average value. If a = 1 (that is, if all instantons and anti-instantons 
form molecules), one recovers the model of Vanderheyden and Jackson. Let us stress that 
there are good arguments to assume that for the physical phase transition a is in fact close 
to 1. We shall discuss some of them in the conclusion. 



In (0), m is the current mass of the quark and t] the source term for the diquark con- 
densate < ipJP^'^i where Pa = ^(^15^2 projects on a colour 3, scalar diquark state (C 
is the charge conjugation matrix), rj has to be taken to zero at the end of the calculations. 
A general hermitian interaction Ti. can be written as an expansion over the sixteen Dirac 
matrices times the colour matrices A'^: 
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K=l a=l 



(3) 



The measure VTi associated with the random matrices A^°- is: 

Vn = {l[I[ VAf:} exp [-AT E PATr[Af:{Af:r]] (4) 

Ka Xk Ka Aft 

with "DAfi^ being the Haar measure. Pk = 1 if A' corresponds to vector or axial- vector 
interaction and = 1/2 if AT is scalar, pseudo-scalar or tensor. The variance T^xa is the 
same for all channels. 

Following |T^,[TB|, one first performs the integration over the random matrix interaction 
and then uses a Hubbard-Stratonovitch transformation to introduce two auxiliary variables 
cr and A, associated respectively with the chiral and diquark condensate. After integrating 
out the fermion fields, one obtains the following partition function: 



a,d) = J dadA exp [-4:Nn{a, A)] 

with the thermodynamic potential Q{a,A): 

Q{a, A) = AA^ + Ba^ -\ [a{Nc - 2) In [((a + m + ^if + d^){{a + m - iif + d^) 

+ (1 - a){N^ - 2) In [(cr + m + ^f{a + m - 
+2a In [((o- + m + + + A^){{a + ni - ^if + d^ + A^) 
+2(1 - a) In [{{a + m + ^if + A^){{a + m - ^if + A^)] } 



(5) 



(6) 



As discussed in ||T2|,|T3|, the couplings B and A are weighted averages of the Fierz coefficients 
obtained by projection of the original interaction on chiral and scalar diquark channels 
respectively. The ratio B/A is the only independent parameter and measures the balance 
between chiral and diquark condensation; varying this ratio allows to explore all the different 
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possible structures of the phase diagram. Imposing the interaction to be hermitian gives 
the upper bound B /A < The absolute magnitudes of A and B play a secondary role. 

They introduce a scale for the condensating fields but don't affect the structure of the phase 
diagram. In the following calculations, we fix A = 1 and vary B. 

To study the various phases for a given value oi B/A one must minimize the potential 
VL with respect to both condensates and solve the resulting gap equations. One can already 
notice that, if a 7^ 1, the second term of this potential is singular at zero density if the chiral 
condensate vanishes, a = (the current quark mass m will be set to zero in the following 
calculations). This is due to the fact that chiral symmetry cannot be completely restored at 
zero density unless all the instantons and anti-instantons form molecules (i.e. a = 1), see 
also |0| . 



III. RESULTS AND DISCUSSION 

The gap equations derived from the potential (|^) admit four kinds of solutions: the 
trivial vacuum where both condensates a and A are zero; the chirally-broken phase where 
cr 7^ and A = 0; the colour-superconducting phase (A 7^ 0,0" = 0); the mixed-phase, 
where both condensates are non-zero. By varying the ratio B /A^ we recovered the various 



scenarios discussed in |]T3[; in particular, the mixed-phase appears only for B/A > 1.05. 



We first present results for the coupling ratio B/A = -{— — - — ) = 0.75 corresponding to 

2 N(, 1 

one-gluon-exchange as well as to an instanton-induced interaction. Fig.|l] shows the phase 
diagram in the ^ — d plane for six different values of the molecule fraction a: 1, 0.99, 0.9, 
0.7, 0.5 and 0.1. Such phase diagrams can be misleading if the actual condensate values are 
too small to be stable against fluctuations. This does, however, not seem to be the case, as 
can be seen from Fig.0 and ^. 

For a = 1, one recovers the results of [|I^] by identifying the parameter d with the first 
Matsubara frequency ttT. We have checked that the restoration of chiral symmetry is first- 
order for small enough d {d < 1.57). The phase diagram changes noticeably for a ^ 1: 
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FIG. 1. Phase diagram for B/A=0.75. The x-axis corresponds to the chemical potential; the 
y-axis to the parameter d. x ^ label the chiral and diquark phase respectively. Dashed and 
continuous lines correspond respectively to first and second-order phase transitions. Note that 
Q CD-phenomenology suggests that at the phase transition a is close to 1. 
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FIG. 2. The strength of the quark condensates a in Fig.! {BjA = 0.75) 
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FIG. 3. The strength of the diquark condensates A in Fig.|^ 
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as anticipated previously, chiral symmetry cannot be restored at zero-density. (This is well 
known, as a 7^ 1 implies the presence of isolated instantons and thus zero-modes.) The chiral 
transition becomes first-order for all values of d, so there is no longer a tricritical point in 
the phase diagram. Moreover the diquark phase appears also at large d and relatively low 
density. If a is further decreased, the diquark phase extends to all values of d and for growing 
values of d the chiral symmetry restoration occurs at higher /x. One notes, however, that 
chiral symmetry restoration always occurs for above a certain critical value, even if the 
fraction of instanton-molecules is small. We conclude that for small /i the properties of the 
chiral phase transition depend crucially on the detailed instanton-anti-instanton dynamics. 
In contrast, the occurrence of a diquark phase seems to be a generic property of the phase 
diagram whatever the fraction of molecules. 

For larger value of the ratio B/A a. mixed phase appears: we show in Fig.^ results for 
the case B/A = 1.4. 



Again we reproduce the results of [|T^] in the case a = 1. The transition from the mixed- 
phase xA to the diquark phase A is first-order for d < 1.1. For a < 1, the changes in the 
phase diagram are similar to the ones observed in Fig.^. The mixed-phase and the diquark 
phase extend over the whole range of d and the transition from the mixed-phase to the 
diquark phase is now first order for all values of d. As in the previous case the diagram is 
stable only for high enough fi. 



For completeness, we have also checked the results obtained in [0 when the current 
quark mass m is non-zero. For a = 1, one recovers their phase diagram, where a first-order 
transition line (corresponding to the first-order chiral transition of the chiral limit) ends in 
a critical point. This critical point however disappears if a 7^ 1 and one finds a first-order 
transition for all d. 

Before to conclude this section, we consider briefly the case Nc = 2. As it has been shown 
in p]^JT3|] , the ratio B/A is necessarily equal to one in that case. At zero chemical potential, 
the potential Q depends on the condensation fields a and A through the combination cr^+ A^. 
As soon as /i 7^ 0, this symmetry is broken and the chiral condensate vanishes; the system 
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FIG. 4. Same as Fig. 1 but for B/A=1.4. labels the mixed-phase. 
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prefers diquark condensation over chiral symmetry breaking. The same conclusions are 
valid for a molecule fraction a 7^ 1: the only difference is that the diquark condensate stays 
non-zero at low /i even for large values of d. 



IV. CONCLUSIONS 



We have investigated the phase diagram of QCD as a function of density and tempera- 



ture within Random Matrix Theory along the lines of ||T2|,|T3[. The aim was to distinguish 
generic and specific properties. We did not treat the temperature dependence explicitely 
but encoded it in the fraction a{fi,T) and strength d{fi,T) of instanton-anti-instanton cor- 
relations. We studied the stability of the phase diagram with respect to variations of a and 
d in comparison with the results of [0. We found that the low-density part of the phase 
diagram is highly dependent on T) and d{fi,T) and thus on the detailed instanton 
dynamics, while, on the other hand, the phase diagram is rather stable for high densities. 
We also found that for finite /i not all instantons and anti-instantons have to combine into 
molecules to allow for a chiral phase transition. These are our main results. Of course. 
Random Matrix Theory only allows us to derive qualitative results. In particular, it can 
not address the question mentioned in the introduction about the dependence of the gap on 
the QCD coupling constant; this dependence is closely related to the behaviour of the gluon 
propagators while our approach amounts to use a contact quark interaction. 

For quantitative studies we need realistic models for the functions T) and d{fi,T) 
i.e. specific input which goes beyond RMT. The fraction of molecules has been calculated 



as a function of temperature at zero density in |21|, where it was shown that the fraction of 



molecules jumps rapidly from a ~ 0.5 at T = 0.7 Tc to 1 at Tc. In P,p2| the molecule density 



was computed as a function of the chemical potential at zero temperature. It was found that 
there is a delicate competition between random instantons engaged in diquarks and molecule 
formation: in fact, the former are dominant and induce the chiral transition. The fraction of 
atomic instantons goes, however, also in this case to zero when chiral symmetry is restored. 
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In our case also, the diquark phase appears in the regions of the phase diagram where the 
effects of molecules are weak enough (meaning either a small fraction a or a relatively low 
value of d) . 

Relating our phase diagrams to the one in the {/i, T) plane is not a trivial question. We 
have shown that the phase diagram is rather unstable at low density, which means that 
its characteristics in the (a*, T) plane should depend on the details of the functions a{ii,T) 
and d(//,T). The results quoted above indicate that the fraction of molecules a changes 
abruptly at the phase transition. It is natural to assume that this will also happen for the 
intermediate case (T and ii non-zero), which suggests a large value of a. Whether this value 
will be exactly equal to one (which, as we have seen, is necessary for a vanishing of the chiral 
condensate) or not is not known. However, if the value of a is large enough (> 0.9), the 
chiral condensate will anyway not be stable against fluctuations. We feel that our present 
understanding of the QCD phase transition is insufficient for reliable quantitative models 
for the functional form of «(//, T) and d(//, T). 

This work was supported by BMBF(GSI) and DFG. We acknowledge helpful discussions 
with M. Alford, A. Jackson, B. Vanderheyden, and H.A. Weidenmiiller. 
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